A generating mechanism for higher order rogue waves 
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We introduce a mechanism for generating higher order rogue waves (HRWs) of the nonlinear 
Schrodinger(NLS) equation: the progressive fusion and fission of n degenerate breathers associated 
with a critical eigenvalue Ao, creates an order n HRW. By adjusting the relative phase of the 
breathers at the interacting area it is possible to obtain different types of HRWs. The value Ao 
is a zero point of the eigenfunction of the Lax pair of the NLS equation and it corresponds to an 
infinitely large period of the breather. By employing this mechanism, we prove two conjectures 
regarding the total number of peaks, as well as a decomposition rule in the circular pattern of an 
order n HRW. 
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Introduction. Rogue waves (RWs) in ocean are catas- 
trophic natural phenomena with a long history and fas- 
cinating mariner stories [1[. Detailed studies of RWs have 
only occurred during the past five decades [2|-[s| • A pro- 
totype one dimensional rogue wave is the so called Pere- 
grine soliton[6]; this soliton exhibits the two remarkable 
characteristics of first order RWs: a) localized behav- 
ior in both space and time, b) existence of one dominant 
peak. RWs have been observed in several fields, including 
optics 0-9L superfiuid Helium [To[, Bose-Einstein con- 
densates p lasmas 12[ [l3|, microwaves jl^, capillary 
phenomena[15j, telecommunication data streams [16| and 
inhomogeneous media [Tt(| . 

A typical modeling equation of RWs in fibre optics is 
the celebrated nonlinear Schrodinger(NLS) equation[l8|. 
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Here q = q(x, t) is a complex smooth function of x and 
t, and the subscripts denote partial derivatives. The 
Peregrine soliton[6|, which is the first order RW[19j of 
the NLS equation, has been observed experimentally in 
fibers [201, m a water tank [21] and in multi-component 
plasma[22j. Recently, a super rogue wave [23], i.e. the 
second order RW, has also been observed in a water 
tank. In addition to NLS equation, the Hirota equa- 
the first type derivative NLS (DNLS) equa- 




the third type DNLS equation[27|, the NLS- 
Maxwell-Bloch equations [28], the discrete NLS equa- 
tion [29], the two-component NLS equations (30l-[32| and 
the Davey- Stewart son equation [33l l34j, also admit RWs. 
These results show that RWs may be generic phenomena 
in nonlinear systems. 

The Peregrine soliton 0, [l9[ of the NLS is expressed 
in terms of a simple rational formula, it corresponds to 
a simple profile, and can be obtained from a breather 
solution via the simple limit of the period of modula- 
tion approaching infinity. However, higher order rogue 
waves (HRWs) [35l - l38| are expressed in terms of compli- 
cated formulae and their profiles exhibit several differ- 
ent interesting patterns [39l443|. These patterns include 
a fundamental pattern consisting of a simple central 
highest peak surrounded by several gradually decreasing 
peaks(see Fig 2a and Fig. 3 in reference [43]), equal- height 
triangular patterns (see Fig 2b of the same reference and 
Figure 2 of reference |40|), and equal- height circular pat- 



tern(see Fig 4 of the same reference). 

Taking into consideration the complexity of the rele- 
vant formulae 44[ 45| , as well as the plethora of the dif- 
ferent possible patterns, it is a challenging problem to 
elucidate the mechanism of HRWs generation. There ex- 
ist two important conjectures regarding HRWs. 

• In the case of a single fundament pattern, an order 
n RW has n(n + 1) — 1 non-uniform peaks [43j; in 
the case when there exist several patterns, an order 
n RW has n(n + l)/2 uniform peaks 36|, |37[. 



• In the case when an order n RW displays a ring 
structure, the outer ring has 2n — 1 uniform peaks, 
and the inner structure is an order n — 2 RW[43j. 

In this letter, we present a generating mechanism for 
HRWs of the NLS equation, and using this mechanism, 
we prove the above two conjectures. Furthermore, we 
discuss several new interesting patterns of HRWs. 
Degenerate n-fold DT and inverse DT. In order 
to study the breather and the RW solutions of NLS, 
we shall use the determinant representation of the Dar- 
boux transformation(DT) introduced in [46l-[48|. Fur- 
thermore, we shall use the notions and the main results 
of these references regarding the n-fold DT (theorem 1 in 
[481 ]) and the related functions (gW y-W^H) generated 
by the n-fold DT (corollary 1 in [48j). In order to sat- 
isfy the reduction requirement = — (r^)*,we choose 
hk = {-fk-i i) T >* = 1,2,--- ,n, where T de- 

notes matrix transposition and asterisk denotes complex 
conjugation. Under this reduction, q^ is a solution of 
the NLS generated by n-fold DT starting with the seed 
solution q. In the following we always use this reduction 
condition. 

Theorem 1 and corollary 1 cited above imply that 
a n-fold DT T n of the NLS, annihilates its indepen- 
dent generating functions, which are the eigenfunctions 
fi(i = 1, 3, 5 • • • , 2n — 1) associated with n distinct eigen- 
values Ai,A3,A5,--- , \2n-i- This means that if we fix 
the given set of eigenvalues, we cannot apply DTs more 
than once. Recall that the formulae for the eigenfunc- 
tions fi(i = 1, 3, 5 • • • , 2n — 1) differ only by the fact that 
they involve different eigenvalues A^. However, in order 
to obtain a HRW for a critical eigenvalue Ao, we must 
apply repeatedly DTs. This difficulty can be overcome 



by noting that the annihilated eigenfunctions can be re- 
created by taking the limit A^ — ^ Ai of the used eigenval- 
ues in DT 0. We set ft = 0(A*) and f\ n] = ^ n] \x=x z - 
We shall use the determinant representation of the n-fold 
DT (see eq.(14) of [Hj]) to illustrate the relevant con- 
struction. It is straightforward to verify that /j 1 ' = 0, 
and hence we cannot apply DT again with the eigen- 

/l 1] (Ai + e) = 

=o )e + 0(e). 

Hence, the limit lim ' f!1J — ;i v 1 ^ 

yields a transformed eigenfunction associated with Ai, 
which can be used to generate a new DT so that we can 
apply this DT with the given eigenvalue Ai for the sec- 
ond time. Similarly, set the second degeneration eigen- 



value Ai. Let A3 = Ai + e, then — J3 
/ W(A 1+ e) = /I 1 '(A 1 ) + ( ^ 11( ^ +e) 

1 f in _ 



value A5 = Ai + e in f t 
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the limit lim e ^o ~^f^ = 



/j 2 ' re-creates a transformed eigen- 



d 2 jf (A 1 + e) , 

function associated with Ai of the 2-fold DT. Note that 
the zero order and the first order terms of e in fj^ yield 
zero contributions. In general, for a n-fold DT, we can 
use the following theorem on (f)^(X) and using the 
degeneration limit A^ — >• Ai by a similar analysis, based 
on the determinant representation given by theorem 1 
and corollary 1 of [48] . The following notations, includ- 
ing matrix elements ((£1)12)17 anc ^ (W^niij, are given in 
0. 

Theorem 1 A n-fold DT with a given eigenvalue Ai is 
realized by the degenerate limit Xi — >• Ai . This degener- 
ated n-fold DT yields the transformed eigenfunction (j)^ 
of X, where 
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as well as a new solution of the NLS equation starting 
with the seed solution q,where 



<7 [n] (a:,i;Ai) =q-2i 
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with 
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6n-i,i(Ai + e) 

£=0 / 2nxl 
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QUi — 1 



6n,i ( Al +C) 
£=0 / 2nxl 



(Q2n)ij(Ai + e) 
£=0 / 2nx2n 



rii = [^jp], [i] denotes the floor function of z, Q 2n is the 
determinant in the numerator of (£1)12 148[ . 
Starting with different seed solutions q, eq.(|3]) yields dif- 
ferent degenerate solitons and breathers. Furthermore, 
by choosing a special eigenvalue Ai = Ao associated with 
0(A O ) = , Eq.© yields an order n RW. 



In order to get an order n — 2 RW from an order n 
RW by a simple limit, it is necessary to use an inverse 
DT. For a general eigenvalue A, the x— part of the 
Lax pair of the NLS admits the solutions 0(A), as well 
as the linearly independent solution ip(\) = (^i,^2) T - 
Furthermore = T n %j) and 0^ = T n (j) are linearly 
independent because T n is a linear transformation of 
^ and <j>. Let g k = (g k i,9k2) T = ^(Afe); then the 
Wronskian determinant W{fi,gi) = fngi2 ~ Mgn of 
fi and gi is a non-zero constant. Using the determi- 
nant representation of the one-fold DT, T(A;/i,/ 2 ), 
generated by fi and f 2l we find the transformed 

[i] _ (A 1 -A 2 )^(/i,^i) 
9\ — 



functions 



(\i-\ 2 )W(f 2 ,g 2 ) 



\W 2 \ 



zero in contrast to /j 1 ' = and f 2 l = 0. Hence, we 
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which are not 



[2] 



can use and g 2 l to generate the second fold DT 
T(X; g^). Using a straightforward calculation with 
the help of Theorem 1 in [48] , it can be shown that the 
two-fold DT is given by T 2 = T(A; g^\g^ ] )T{X- f u f 2 ) = 
(A — Ai)(A — A 2 )/. Here we just present the cal- 
culation of the element (T2)n. First note that 
\W 4 (g 1: g 2 Jij2)\ = -(A 2 - X 1 ) 2 W(f 1 ,g 1 )W(f 2 ,g 2 ); 
(T 2 )n = -(X-X 1 )(X-X 2 )(X 2 -X 1 ) 2 W(f lj g 1 )W(f 2j g 2 ). 

(r 2 )n 



Hence, (T 2 )n 



(A - Ai)(A - A 2 ). 



W A {g 1 ,g 2 J 1 J 2 ) 
Thus T(A;#[ 1] ,^ 1] ) is the inverse DT of T(A;/i,/ 2 ). 
In general, for an (n-2)-fold DT T n _2 generated 
by /i,/ 2 , • • • ,/ 2 n-5,/2n-4, we can find one-fold in- 
verse DT as follows(note that g 2 n-3 = ^n-iQin- 



-3 or 



[ n -2] _ t> 
92n-2 — 1 n-2Q2n-\ 

independent): 
Theorem 2 Let 



and /2n-3° r f [ 2 U n- 



2 are linearly 



(2) T(X;ftzlft- 2 ] ) 



the 

1 2n-3 •> J 2n-2 ) a ^ er 

and q [ Z~-l = T(X;f 2n 

KJn-2] 
-2 )92n-2 ' 

T(\- J n ~^ J n ~^ 
1 \ A i92n-3i92n-2 

DT. 



(n 
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T(X'J 2n 



[n-2] ,[n-2]^ 
3w2n- 



is the inverse of the (n 



- l)-th fold DT be 
(n-2)-fold DT T n _ 2? 

2] [n 
3 J ^2n-2 

n-th fold DT 
D-th fold 



n-2] ,[n-2k [n-2] ^[n-1] 
-3 'J2n-2 Ji/2n 

Then 



In other words, the n-th fold DT T(A; ^ 2 n-3 j ^2n-2 ) 



maps <7 L 



to g 



[n-2] 



This gives an important connec- 



tion between RWs of order (n-1) and order (n-2). 
Higher order breathers and rogue waves The 

first order breather of the NLS equation is a periodic 
traveling wave. This solution, via the limit of the period 
approaching infinity, gives the first order RW p, fl9| . 
However, it is still not clear how to generate HRWs 
from multi-breathers, even for second order RWs (4lj . 
Moreover, the collision of three breathers [49] does not 
provide a satisfactory explanation for the appearance of 
different patterns of order 3 RWs. 

On the (x, t) plane, because of the conservation of the 
number of the breathers, there exist n separate peaks in 
each row before and after the interaction of n breathers. 
The interaction area is localized near the origin of the 
plane between the two closest rows( or periods) pos- 
sessing n separate peaks. When the breathers are in 
the interaction area, their peaks get closer. Based on 
the detailed investigation of the interaction of breathers, 
we claim the following mechanism for the generation of 



HRWs: the progressive fusion and fission of n degenerate 
breathers associated with a critical eigenvalue Ao , creates 
an order n HRW. Furthermore, by adjusting the relative 
phase of the breathers at the interacting area, it is pos- 
sible to obtain different patterns of HRWs. Here Ao is 
a zero point of the eigenfunction 0(A), i.e. 0(Ao) = 0, 
which corresponds to the limit period of breathers be- 
comes infinitely large. The relative phase can be adjusted 
via the tuning of the parameters Si in the eigenfunctions 

In this Letter, we shall take a periodic seed q = ce lp 
with p = ax + (2c 2 — a 2 )t. The corresponding eigenfunc- 
tions 0(A) = (0i,02) T is given by 

,^ w /ceH«+^)) +i (| + Cl (A)+A)e-H-f+^))\ 
^ Aj "V^- i( « +d(A)) + i(f + c 1 (A)+A)e i (-«+^)) )• W 

Here Ci (A) = ^c 2 + (A + a/2) 2 , d(A) = c ± {X){x + (2 A - 
a)t + s + $), $ = X^=i 5/ee 2/c [40], n denotes the num- 
ber of the steps of the multi-fold DT, Ao = —a/2 + ic 
is a zero point of the eigenfunction 0(A), e denotes a 
small parameter when we consider the degeneration of 
the eigenvalues, i.e., A = Ao + e, are complex constants. 
The functions fa = 0(A^) have the same form except for 
different values of the eigenvalue which is necessary to 
generate HRWs via the process of the eigenvalue degen- 
eration Xi \-> Ai (see Theorem 1). In the following exam- 
ples we set a = 0. Also, in order to adjust the relative 
phase of the breathers at the interaction area according 
to Theorem 1 and corollary 1 of reference (isj . we set 
Si = 0(i > 1), but so has different values in different fa. 
There exist three types of relative phases of n breathers 
in the interaction area: synchronous, anti-synchronous 
and quasi synchronous. 

In the interaction area of n synchronous breathers, 
there exist progressively increasing fusion via n — 1 steps 
from the n lower peaks to the central maximum peak, and 
then progressively decreasing fission via n — 1 steps from 
the central maximum peak to the n lower peaks. Here, 
each step of fusion annihilates one peak and hence the 
height of peaks increases; similarly, each step of fission 
creates one new peak and hence the height of peaks de- 
creases. These peaks are arranged as two triangles with 
one joint vertex along their perpendicular bisector. Thus, 
the total number of the non-uniform peaks in the inter- 
action area is n{n + 1) — 1. It is interesting to note that 
the outermost row of the interaction area has n lower 
peaks, which are closed to each other. Hence, peaks are 
much lower than the ones in the nearest row of the non- 
interaction area. This phenomena provides evidence for 
the strong interaction of the breathers. When the eigen- 
value used in the breathers approaches the critical value 
Ao,i.e. Xi ^ Ao, the period of all breathers goes to infin- 
ity simultaneously, so that the only profile in the interac- 
tion area survives, which gives the fundamental pattern 
of a HRW. Therefore, this pattern of an order n HRW, 
has n(n + 1) — 1 non-uniform peaks. The central profile 
of the three breathers in Figure 1 is very similar with 
the fundamental pattern of an order 3 RW(see Fig3(a) 
of [43[) of the NLS equation. The three breathers are 
plotted according to Theorem 1 and corollary 1 of (48| 
with a = 0.01,c = 0.5, s = 0,and Ai = -0.2 + 0.54i in 
/1, A 3 = 0.1 + 0.55i in fa and A 5 = 0.03 + 0.56i in fa. 



The interaction of n anti-synchronous breathers is 
simpler, although there exists also the fusion or fission 
of peaks. In the interaction area, peaks are closer to 
each other. By suitably adjusting the relative phases 
of breathers, n synchronous breathers become n anti- 
synchronous breathers, and the corresponding peaks in 
the triangle disappear, so that only peaks in one trian- 
gle survive. Specifically, by suitably changing the rela- 
tive phase, we observe the disappearance of the lowest 
peak in the outermost row of the interaction area, then 
the two nearest peaks disappear, and so on. This chain 
reaction continues until the coalescence of the two trian- 
gles. The collapse of this triangle is stimulated by the 
loss of the nearest-neighbour interactions. Thus, there 
are n(n + l)/2 peaks in the interaction area, which are 
allocated on the remaining triangle. Furthermore set 
Ai >• Ao simultaneously; the profile in the interaction 
area of the order n breather yields a triangular pattern 
of a HRW. Therefore, there are n(n + l)/2 equal-height 
peaks in the triangular pattern of an order n HRW. A 
triangular structure of the 3 anti-synchronous breathers 
is plotted in Figure 2 by using Theorem 1 of [48] with 
a = 0.01,c = 0.5, and Ai = 0.05+0.531i and so = 16 in fa; 
A 3 = 0.55i and s = -20i in fa; A 5 = -0.05 + 0.551i and 
so = — 16 in fa. By suitably choosing different values of 
the parameters in n anti-synchronous breathers, so that 
the relative positions of peaks are changed but the total 
number of peaks is preserved, we obtain a ring structure 
associated with the n(n + l)/2 peaks in the interaction 
area, which gives rise to the circular pattern of an order n 
RW in the above limit possessing n(n + l)/2 peaks. Fig- 
ure 3 confirms the ring structure of 3 anti-synchronous 
breathers with a = 0.01,c = 0.5, and Ai = 0.05 + 0.54i 
and so — 1 + i in /1, A3 = 0.55i and so = in fa, 
A5 = —0.05 + 0.56i and so — 1 + i in fa. There exist 
many other patterns appearing at the interaction area 
of n anti- synchronous breathers, which gives rise to sev- 
eral types of HRWs such as a pentagram (Figure 4 for 
an order 5 RW given by Theorem 1 with Ao = ic, a = 
0,c = Mo = 0,si = 0,s 2 = 1000000,s 3 = 10, s 4 = 0) 
and a triangle in circle (Figure 5 for an order 5 RW given 
by Theorem 1 with Ao = ic, a = 0,c = Mo = 0, si = 
15, s 2 = 0, s 3 = 0, s 4 = 10000000). 

By suitably choosing the values of the parameters in 
fa, we observe n quasi synchronous breathers from theo- 
rem 1 and corollary 1 in [48|. There exist many different 
patterns in interaction area, which implies many inter- 
esting types of HRWs when A^ \-> Ao- For example, one 
can find the following interesting decomposition of an or- 
der n HRW: an order n — 2 RW surrounded by 2n — 1 
peaks [IK . The first nontrivial example of this decompo- 
sition is given by the interaction of 4 quasi synchronous 
breathers. Unfortunately, we are not able to plot the 
profile in the interaction area in this case, due to the 
complexity of order 4 breathers. However, we find thor- 
ough decomposition of the circular pattern, an order 5 
RW in Figure 6 with Ao = ic, a = 0, c = Mo = si = 
S2 = 0, s 3 = 80000, s A = 20000000 and an order 6 RW 
in Figure 7 with Ao = ic, a = 0,c = 1,sq = s\ — ,$2 = 
0, s 3 = 30000, 54 = 0, s 5 = 100000000. 

The inverse DT provides a technique enabling us to 
prove the above interesting decomposition rule of an or- 



4: 



der n RW. For simplicity we set a = in the seed so- 
lution, and then we set Ao = ic. Expanding 0(Ao + e) 
with respect to e, we find that the coefficient of the 
first order of e is given by fo = (e lc t (2\cx — Ac 2 t + 
2\csq + i), — e~ lc t (2icx — 4c 2 t + 2ics — i)) T , which is 
an eigenfunction associated with q and Aq. There exist 



another eigenfunction go = (e 



-\c 2 t\T 



Y of A . Note 



that /o and #o are two linearly independent eigenfunc- 
tions of Ao- According to Theorem 1, an order n RW 
is generated by a degenerate n-fold DT with the critical 
eigenvalue Ao from the periodic seed q = ce 2lc . Let 
T n _i be a (n-l)-fold degenerate DT with Ao, so that 
g n-1 ^ and f^ 1-1 ^ are given by eq.@, and are linearly 
independent. By a tedious asymptotic analysis, we find: 

/d n_1 ^ = fo + Sn-i^o 1-1 '' wnere fo is a smooth bounded 
function. According to Theorem 2, the n-th fold DT de- 
fined by T(A; g [ n ~ 1] ) is the inverse of the (n - l)-th DT 
defined by T(A;/q U_ ^). Thus, by the limit s n -i \-> oo, 
the n-th fold DT T(X;fl n ~ 1] ) = T(A;# n-1] ), gives an 
inverse transform of the (n — l)-th fold DT. Therefore, 
under this limit, an order n RW q^ n \ reduces to an order 
n — 2 RW, q\- n ~ 2 \ By taking s n -i to be large (but finite), 
an order n RW is decomposed into an order n — 2 RW 
and 2n — 1 peaks located at an outer circle such that the 
total number n(n + 1)/2 of peaks can be realized either in 
a triangular pattern or in a circular pattern. The inner 
order n — 2 RW, can take any of these formulae by choos- 
ing Si(i = 0, 1, • • • , n — 3). This decomposition rule of 
HRWs was conjectured by Akhmediev' s group [43|. Fig- 
ure 5 and Figure 6 show different patterns of the inner 
lower order RW decomposed from an order 5 RW. The 
fundamental pattern of inner order 3 RW reduced from 
an order 5 RW is shown in Fig. 4(c) of [43[- According 
to this decomposition rule, Figures 6 and 7 provide the 
first two non-trivial examples of a complete decomposi- 
tion associated with three levels. 

Conclusion A new generating mechanism for the cre- 
ation of HRWs for the NLS is introduced based on the 
fusion and fission of n degenerate breathers at the critical 
value Ao- The different types of HRWs are determined 
by the relative phase of the breathers at the interaction 
area, and this can be adjusted by the free parameters 
Si appearing in the eigenfunctions fa. Furthermore, two 
conjectures related to HRWs are verified using the deter- 
minant representation of a degenerate DT and its inverse. 
Our results elucidate the mechanism responsible for the 
generation of various types HRWs, hence they should be 
useful both for the generation, as well as for the preven- 
tion of HRWs in application. 
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Fig. 1: (Color online) The fusion and fission of three 
synchronous breathers on (x,t) plane. The right panel is 
a local profile in the interaction area of the left panel 




Fig. 2: (Color online) The fusion and fission of three 
anti-synchronous breathers on (x,t) plane. The right 
panel is a local triangle pattern in the interaction area 
of the left panel. 




Fig. 3: (Color online) The fusion and fission of three 
anti-synchronous breathers on (x,t) plane. The right 
panel is a local circular pattern in the interaction area 
of the left panel. 




Fig. 4: (Color online) The pentagram pattern of an 
order 5 RW. There are three concentric circles, and each 
of them have five peaks. 
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Fig. 5: (Color online) A triangle pattern in a circle for 
an order 5 RW. The right panel is a local central profile 
of the left panel. 





Fig. 6: (Color online) Decomposition of an order 5 RW. 
The right panel is a local central profile of the left panel. 
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Fig. 7: (Color online) Decomposition of an order 6 RW. 
The right panel is a local central profile of the left panel. 



